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ABSTRACT 

We present a combined analysis of all publicly available, visible HST observations 
of transits of the planet HD 209458b. We derive the times of transit, planet radius, 
inclination, period, and ephemeris. The transit times are then used to constrain the 
existence of secondary planets in the system. We show that planets near an Earth 
mass can be ruled out in low-order mean-motion resonance, while planets less than 
an Earth mass are ruled out in interior, 2:1 resonance. We also present a combined 
analysis of the transit times and 68 high precision radial velocity measurements of the 
system. These results are compared to theoretical predictions for the constraints that 
can be placed on secondary planets. 
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1 INTRODUCTION 

The discovery of 51 Peg B jMavor fc Quelozlll995h fulfilled a many-dec ade quest to discover extrasolar planets arou nd main 
sequence stars. The discovery of the first transiting planet HD 209458b jcharbonneau et aljEnol : iHenrv et alJl200tj> allowed 
the confirmation that radial velocity variations were indeed due to planets and not some other period stellar phenomenon. 
Currently ten transiting planets have be en found, allowing a measurement of the mass-radius relation for short period extra- 
solar planets Jcharbonneau et alJl20o3) . Besides confirming that radial velocity variations are due to real objects with the 
sizes and masses of planet s, transiting planets have allowed a host of other questions to be addressed, for example, the plane t 
atmospheric composition Jcharbonneau et alJl20o3^. constraints on the presence of planetar y satellites terown et al.ll200lth 
detection of thermal emission from the planets ijPe.rm ng et alJ200i4 Icharbonneau et al]l2005h . and the size of the planet core 
as a function of stellar metallicity l|Ouillot et alJl200fiF) . 

In this paper we provide another application of transiting planetary systems: we place a constraint on th e presence of 
terrestrial-mass planet companions in mean- motion resonance with HD 209458b. As shown in lAgol et alJ J2005h . a planet in 
resonance induces a libration of the transiting planet with an amplitude ~ 0.2j _1 (M2/Mi)Pi, where Mi, Pi are the mass and 
period of the known transiting planet and M2 is the mass of the companion planet in a j:j + 1 resonance. For HD 209458b, 
this translates into a variation in the times of transit with a 3 minute amplitude and few month period for M2 = m®. Since 
this technique relies on a precise measurement of the times of transit, we use the highest precision photometry av ailable for 
thirteen transits measured by the Hubble Space Telescope jBrown et al.ll200ll:ISchultz et aljEool iKnutson et alJboOri h We 
carry out an analysis of the transit times with a careful quantifying of the timing errors We use these results to constrain 
the transiting planet parameters (©• We then use these derived times of transit to place a constraint on a secondary coplanar 
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Table 1. Overview of the HST observations of the HD 209458 system including the instrument, wavelength range, cadence, and date of 
observation. 



Transit 


Transit number 


Orbits 


Instrument 


Filter/Grating 


wavelength range 


exposure/readout time 


date 


1 





5 


STIS 


G750M 


581.3-638.2 nm 


60/20 sec 


UT 2000 April 25 


2 


1 


5 


STIS 


G750M 


581.3-638.2 nm 


60/20 sec 


UT 2000 April 28-29 


3 


3 


5 


STIS 


G750M 


581.3-638.2 nm 


60/20 sec 


UT 2000 May 5-6 


1 


5 


5 


STIS 


G750M 


581.3-638.2 nm 


60/20 sec 


UT 2000 May 12-13 


5 


117 


3 


FGS 


F550W 


510-587.5 nm 


0.025/0 sec 


UT 2001 Jun 11 


6 


143 


2 


FGS 


F550W 


510-587.5 nm 


0.025/0 sec 


UT 2001 Sep 11 


7 


1(50 


2 


FGS 


F550W 


510-587.5 nm 


0.025/0 sec 


UT 2001 Nov 10 


8 


179 


2 


FGS 


F550W 


510-587.5 nm 


0.025/0 sec 


UT 2002 Jan 16 


9 


252 


2 


FGS 


F550W 


510-587.5 nm 


0.025/0 sec 


UT 2002 Sep 30 


10 


313 


5 


STIS 


G430L 


290-570 nm 


22/20 sec 


UT 2003 May 3 


11 


321 


5 


STIS 


G750L 


524-1027 nm 


19/20 sec 


UT 2003 May 31 


12 


328 


5 


STIS 


G430L 


290-570 nm 


22/20 sec 


UT 2003 Jun 25 


13 


331 


5 


STIS 


G750L 


524-1027 nm 


19/20 sec 


UT 2003 Jul 5-6 



planet fSHjl. We finish up with a comparison of these results with those derived for the TrES-1 transiting planet system ( H4.H 
and a comparison with the constraints from radial- velocity measurements f >!4.2ll . 



2 DATA ANALYSIS 

At least ten science observing programs with the Hubble Space Telescope have proposed to observe transits of the planet 
in the HD 209458 system. The five ultraviolet observations do not have enough photons to measure the transit lightcurve 
precisely, while four programs have made observations of the planetary transit in the optical, of which three are now publicly 
available: 8789, 9171 and 9447 (the other program, 10145, has not been fully made public yet). Two programs, Brown 8789 
and Charbonneau 9447, utilized the STIS spectrograph as a CCD photometer. By spreading the light from the star with a 
grism, in an effort to capture as many ph otons as possible, their relative photometry is likely the most precise ever obtained 
jBrown et ailbOQli IrCnutson et~ai] 2006). The third program utilized the Fine Guidance Sensors to carry out rapid high 
precision photometry JSchultz et alj|2004h . These programs give a total of thirteen transits observed with HST. The primary 
motivation for these observations was to constrain the planetary parameters, to detect absorption features in the planet's 
atmosphere, and t o search for satellites of t he planet. However, these data are also useful for timing each transit to measure a 
precise ephemeris JWittenmver et al1l2005l) and to look for deviations from exact periodicity. A summary of the observations 
is given in table Q 



2.1 Reduction of pipeline data 

We re-reduced these da ta from the p i peline cali brated produc t s tha t are available in the archive. We followed the reduction 
procedures described in lBrown et alJ feOQll) and lSchultz et all i2004l) . but summarize the most important steps for each data 
set here. For the STIS data we carried out the minimum reduction necessary to obtain precise photometry. We subtracted 
cosmic rays using a 5 sigma rejection of the time series for each pixel within an HST orbit, we binned the photons within 
8 pixels of the spectrum peak (17 pixels wide) to derive the total counts for each frame, and assigned the midpoint of each 
exposure as the time for each frame. We discarded the first frame from each orbit. For the first transit observed in program 
87 89 we only used the data for which the spectrum centroid was greater than 4 pixels from the edge of the CCD (see discussion 
m iBrown et aTlboOll) and we did not utilize the first and last columns of the CCD which show larger errors. The 1-sigma 
errors on the flux are taken from the pipeline calibrated errors which are summed in quadrature. For the FGS data we took 
the sum of the four photometers, discarded data which deviated more than 3 sigma from a fourth order fit to the time series 
for each orbit, and binned the data to 80 second bins (10 second bins gave the same results) with the time stamp at the center 
of each bin. With these time series we then fit a model for the transiting planet and the flux sensitivity variations of Hubble, 
which we describe next. 
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2.2 Photometric Error bars 

The errors computed from the pipeline include counting noise and, for the CCD data, read noise. We examined the scatter 
in the data compared to the computed error bars and found that the scatter was larger than the pipeline errors for the FGS 
data and second set of STIS data. This could not be attributed solely to cosmic rays and we were not able to identify the 
additional error source. Consequently, we estimated the errors for each transit individually by analyzing the out-of-transit 
data. 

To account for the fluctuations in sensitivity of HST when we determine the size of the error bars, we fit a model to the 
out-of-transit data for each transit with two components: i) an nth order polynomial fit as a function of the orbital phase 
of HST ii) a linear fit as a function of time. We measured the standard deviation of the residuals of this fit in units of the 
statistical error of each point as a function of n, and measured the median of this standard deviation for all thirteen transits as 
a function of n. We found that between n — 5 and n = 6 the median of the standard deviation ratio was constant, indicating 
that including a 6th order term did not significantly improve the fit to the fluctuations in sensitivity as a function of the HST 
orbital phase. The residuals did not correlate with the HST orbital phase, so we then assigned the greater of the standard 
deviation of the residuals for n = 5 or the statistical error. The first five transits (Table had statistical errors nearly equal 
to the residuals, while the error bars of the remaining transits were increased by 20-50%. These changes to the error bars were 
used throughout the remainder of the analysis. 



2.3 Planet orbit and transit model 

We model the lig htcurves by treat i ng the planet and s tar as perfectly spherical, and by assuming that the planet is on a 
circular orbit fLaughlin et al . 2005; D eming et al.ll2005f) . Thus, there are only three free parameters to describe the planet's 
orbit: the ratio of the planet to star radius, R p /R 3 ; the ratio of the relative velocity of the planet and star to the stellar radius, 
v/R 3 ; and the inclination of the planet's orbit, i. In addition, we allowed the time of each transit to be a free parameter, 
giving thirteen transit times. From these is derived the period of the planet, P, which then gives the ratio of the semi-major 
axis to the stellar radius, a/R 3 = (v/R 3 )P/(2n), necessary for computing the relative position of the planet and star. In 
this approach the actual values of the planet and stellar radii cannot be determined, only their relati ve values (to determine 
the absolute values requires fitting the radial velocities and using a stellar mass-radius relation as in lCodv fc SasselovlbooJ : 
iKnutson et"afll2006h . 



2.4 Limb darkening 

If the lightcurves for each transit were only affected by Poisson noise and were continuously sampled, then assumptions about 
the shape of the transit lightcurve would be unnecessary: the mid-point of transit could be found by simply assuming the 
transit lightcurve was symmetric about mid-transit. However, the HST data are affected by sensitivity variations due to 
temperature changes in the telescope and by observing gaps, both due to Earth occultation. This fact forces us to choose a 
model for the transit lightcurve in order to solve for the transit time, i.e. time at transit midpoint. The primary uncertainty 
in the modeling is the shape of the stellar limb darkening. To test the robustness of the transit times to assumptions about 
limb-darkening, we modeled the limb-darkening of the star with two approaches: stellar atmosphere predictions and quadratic 
limb-darkening. 

I n the first case we used the computed limb-darkening from Phoenix NextGen stellar atmosphere models ijciaret fc Haxiscrfij 
2003). These models include an accurate, high-resolution treatment of the brightness at the limb using spherically-symmetric 
stell ar models. We used an atmosphere of 6100 K and a surface gravity of log(g) = 4.5, appropriate for HD 209458a 
JValenti fc Fischeill2005h . We convolved the wavelength-dependent limb-darkening with the number of photons detected at 
each wavelength (for the FGS we convolved the FGS sensitivity with a stellar atmosphere to estimate the relative number 
of photons detected as a function of wavelength). The resulting limb darkening is shown in Figure Q There is only a grad- 
ual change in the limb darkening between the four different instrument/filter combinations, but given the timing precision 
required, these differences are significant. We then used a cubic spline interpolation of the limb darkening to integrate over 
the portion of the stellar disk occulted by the planet. The result is a more gradual ingress and egress than that predicted by 
analytic limb-darkening models which have a sharp edge. 

In the second case we assumed that the limb-darkening was a quadratic function in the cosine of the polar angle measured 
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Figure 1. Limb darkening for the four instrumental configurations - the vertical axis is the surface brightness of the star, /, as a function 
of fi = cos0, where 8 is the polar angle measured from the line of sight to the star. Solid line is for STIS G750M, dotted line is for FGS 
F550W, dashed line is for STIS G430L, while dash-dot line is for STIS 750L. 



from the line of sight. We utilized the routines from iMandel fc Agoll |2002) for computation of the lightcurves (this led to 
eight free parameters for the four filter/instruments combinations used in the data). We primarily used this case to check 
that the timing results are insensitive to the accuracy of the limb-darkening predicted from the atmosphere model, discussed 
more below. 



2.5 HST flux sensitivity variations 

Thermal fluctuations during the orbit of HST lead to changes in the sensitivity of each instrument at the 0.5% level. Fortunately 
these variations appear to be fairly reproducible during successive orbi ts (in addition to a s mall er linear drift betwee n successive 
orbits), but they are not reproducible on longer timescales. Following iBrown et alJ ll200lh and lSchultz et al l J2004l we phased 
the observations of each transit to the HST orbital period (ph)- We used two models for the variation of the flux sensitivity 
with orbital phase and two models for the "secular" variation with time over each transit, combined together giving a 
total of four models for the background variations. The orbital phase variation we modeled as i) a fifth-order polynomial, 
y^._ d(j)^ +1 \ where ^ <f> < 2n is the orbital phase of Hubble, or ii) as a linear and two harmonic functions of orbital phase, 
Co0 + ci sin(ei) + C2 cos(</>) + C3 sin(2</>) + C4 cos(2</>). The secular variation with time we modeled either as a linear function, 
C5ftran(t) + c§t, where ftran is the transit lightcurve model normalized to unity outside transit, or a constant coefficient for 
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each HST orbit, Cj+sJ %ran,i(t) where i labels each HST orbit. The lightcurve model depends on sixteen non-linear parameters 
(the 3 planetary orbital parameters and 13 transit times). For the models of the HST flux sensitivity variations there are 
either 7 x 13 = 91 or 5 x 13 + 51 = 116 linear parameters. Coincidentally, the physically interesting parameters are non-linear, 
while the uninteresting parameters are linear. 

To most efficiently find the best-fit for all 107 or 132 parameters, we used the following procedure: 

1) We ca rried out a minimization over the 16 non-linear parameters using the MINPACK non-linear least squares routine 
jCawerfll984h . 

2) The model lightcurve for each transit we computed by carrying out a linear least-squares fit for the 91 or 116 linear 
parameters. 

This procedure converged rapidly to the best \ 2 ■ We mapped the \ 2 as a function of each non-linear parameter (while 
marginalizing over the other parameters), finding that the \ 2 has a quadratic shape near the minimum, indicating that there 
is a unique solution. 



2.6 Timing Error computation 

The best-fit models have reduced chi-square in the range \ 2 = 1-07 — 1.25 for the four different flux sensitivity models. The 
reduced chi-squares for the FGS data (which contain fewer data points) were as high as 2.6 (for the 5th transit) indicating 
that the models were a poor fit to the data. We examined the residuals of the FGS data and found that some residuals within 
the transits were correlated. Thus our best fitting models were not a sufficiently good description of the data (either due to 
errors in the limb darkening or inaccuracies of the modeling of the HST flux sensitivity variations). Therefore, to obtain a 
better estimate of the error bars, taking into account these systematic effects, we used a Monte Carlo bootstrap simulation 
of the errors in which we shifted the residuals by a random number of points for each transit so that the correlations were 
maintained. We also reversed the residuals for each orbit and repeated the shift. We followed this procedure 200 times, added 
the shifted residuals to the best-fit model, and then re-fit using the same procedure that was applied to the original data 
described in the previous section. Our expectation is that this procedure gives an estimate of the errors due to inaccuracies 
in our assumed flux sensitivity models. We found that the resulting errors on non-linear model parameters, such as transit 
times, were larger by a factor of a few than if we simply used Poisson errors for the bootstrap or resampled the residuals 
randomly. 

Although this process gives an indication of the size of the systematic errors, we also assessed the systematic errors by 
comparing the results from the four different HST flux sensitivity variation models. The differences in \ 2 between the four 
different models was small - all had reduced \ 2 near unity, so it is difficult to favor one particular sensitivity variation model. 
However, the differences between the derived transit times for each model was larger than the size of the error bars for a 
given model. We believe that this discrepancy is probably due to imperfections in all of the four flux sensitivity models, so 
we have taken the mean and standard deviations of all eight hundred Monte-Carlo simulations as an estimate of the transit 
times and their errors. Finally, to check the robustness of the stellar atmosphere limb-darkening we have used the quadratic 
limb-darkening law with the polynomial function of phase and linear secular term to fit the transit times. We find that these 
times agree within the errors with the times derived using the stellar atmosphere model and the same flux sensitivity model, 
so we conclude that the stellar atmosphere limb-darkening model provides an accurate estimate of the transit times of the 
data. 

To these measured transit times we apply a correction for the motion of the Earth about the barycenter of the solar 
system (a heliocentric correction is insufficient as it differs from barycentric by up to 5 seconds). The motion of the Hubble 
Space Telescope about the Earth contributes a negligible timing error. The resulting transit times are reported in Table [5] 
We find from the Monte Carlo simulations that there is no significant correlation between the individual transit times so 
we report the stan dard deviation of ea ch time without the full covariance matrix. Our fitted transit times agree well with 
those obtained by JKnutson et alJfeoOfih . most being within 1/ 2 a. The largest discre pancies are for transits 12 and 13 which 
disagree by 1.25 a and 0.75 a respectively. In the analysis bv lKnutson et all <l2006l) it is these same transits which deviate 
most significantly from a constant period. 
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Table 2. Transit time and uncertainty in days and in seconds for the 13 observed transits. 



Transit 


Transit time 


Error 


Error 




(- 2450000 BJD) 


(days) 


(seconds) 


1 


1659.93678 


0.00015 


12 


2 


1663.46150 


0.00025 


21 


3 


1670.51102 


0.00013 


11 


4 


1677.56044 


0.00037 


31 


5 


2072.33284 


0.00049 


42 


6 


2163.97574 


0.00028 


24 


7 


2223.89685 


0.00058 


50 


8 


2290.86741 


0.00051 


11 


9 


2548.17311 


0.00029 


24 


10 


2763.18306 


0.00018 


15 


11 


2791.38087 


0.00013 


10 


12 


2816.05431 


0.00021 


18 


13 


2826.62867 


0.00017 


14 
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Figure 2. Differences of transit times from uniform period versus number of transit (each point is shifted horizontally slightly for clarity). 
Green/triangle- polynomial and linear; yellow/circle - harmonic and linear; dark blue/diamond - polynomial and flux offset for each orbit; 
light blue/square - harmonic and flux offset for each orbit. Red/cross points are mean and standard deviation of all four flux sensitivity 
models. 
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3 RESULTS 

The results of the best fit models are shown in Figure 
3.1 Constraints on radius, inclination 

Since the data consist of a time series, the best constrained parameters are the period, P, and the duration of each transit, 
T f . This leads to a degeneracy between the physical parameters: inclination and velocity. The relation between velocity and 
inclination is given by 

/ T 2 p 2 \ -1/2 

As long as v and i obey this relation, P and T t are unchanged. Thus, for a very small uncertainty in T t and P, there can 
be a larger uncertainty in i and v/R*. The inclination is constrained by the shape of the limb-darkening and the duration of 
ingress and egress, which partially (but not completely) break the velocity-inclination degeneracy. In Figure 2] we show the 
results of our fitting for v = v/R* and i. We plot the results of Monte Carlo simulations of each of the four models discussed 
above. The correlation between velocity and inclination follows equation Q as expected. 

We derive a radius of R p /R* — 0.1201 ± 0.0006, consistent at la with the radius derived in lMandel fc Agoll ll200l) (our 
derived error is larger since we have taken into account correlated photometric errors using the Monte Carlo analysis). We 
derive a velocity of v/R* = 15.835 ±0.065 days - , which combined with the period gives the ratio of the size of the semi- major 
axis to the size of the star, a/R, = 8.883 ± .036. We derive an inclination of i = 1.5157 ± 0.00095 or i = 86.845° ± 0.055° 
which is consistent with that reported in lBrown et alJ <l200ll) . although our derived error bar is smaller due to our using a 
fixed limb-darkening law. Since one can only derive dimensionless parameters or parameters with dimensions of time or flux 
directly from the lightcurve, we cannot measure the absolute size of the star, orbit, or planet. 



3.2 Ephemeris 

The accurate and realistic estimation of errors on the transit times allows for an accurate estimate of the ephemeris of the 
transiting planet. We find values of: T = 2451659.93677 ± 0.00009 BJD (error is 8 seconds) and P = 3.5247484 ± 0.000000 4 
days (error is 0.03 seconds). Our derived ephemeris is in good agreement with the ephemeris derived bv lKnutson et all J2006I) . 
differing by less than 1/2 a for both P and Tc- We note, however, that the small uncertainty in the period should be considered 
the uncertainty in the mean period that spans the modest time scale over which the data were obtained. Current constraints 
on secondary planets in this system cannot rule out bodies that could change either the mean period or the instantaneous 
period of HD209458b by more than the quoted error. These changes could happen over time scales longer than the span of the 
observations, such as a very small companion in a mean-motion resonance, or over shorter time scales, such as a companion 
on a large, highly eccentric orbit. 

The period and zero-point of the ephemeris are anticorrelated with a correlation coefficient of -0.8; this is due to the fact 
that an increased zero point requires the period to shrink to go through the later data points. The anti-correlation decreases 
to -0.2 if the zeroth eclipse is chosen to be halfway between the first and last eclipses. Given the estimated statistical plus 
systematic errors, the chi-square is y/ 2 = 6.0 for eleven degrees of freedom (thirteen transits minus two fitting parameters). 
This indicates that the transit times are consi stent with bein g unifo rm and that there may not be a second planet detectable 
with the current HST data, as concluded bv lKnutson et all J2006). In the next section we use these data to constrain the 
masses of additional planets in the system. 



4 CONSTRAINTS ON SECONDARY PLANETS 

Following the procedure outlined in Isteffen fc Ago! l|2005l) we use the measured transit times to constrain the presence of 
secondary planets in the HD209458 system. In brief, for a given value of the semimajor axis and eccentricity of a putative 
second coplanar planet we determine the maximum mass that a companion could have without causing a significant deviation 
(3-cr) from the data. Figures [S] and |S| show the maximum mass that a perturbing planet can have in this system for an interior 
perturbing planet and an exterior perturbing planet respectively. 
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Figure 3. Best-fit models overplotted with data (some data -0.15 days before mid-transit is off the plot). The data points are plotted 
with error bars, while the best-fit model is plotted as a red line. Each transit has been normalized to one out of transit and shifted by 
0.05; the bottom is transit 1 and the top is transit 13. 



HD 209458b transits 9 



t — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — r 




86.6 L_i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i 

15.7 15.8 15.9 16.0 16.1 

v/R < (day"') 



Figure 4. Velocity versus inclination - Monte Carlo (green/triangles - polynomial and linear; yellow/circles - harmonic and linear; dark 
blue/diamonds - polynomial and flux offset for each HST orbit; light blue/squares - harmonic and flux offset for each HST orbit) and 
analytic expectation (black line). 



We see from these figures that near several mean-motion resonances, a secondary planet must be comparable to or less 
than the mass of the Earth for various values of the eccentricity. This is particularly manifest in the 2:1 and 3:2 mean 
motion resonances where this technique is most sensitive. Here, a secondary planet can be no larger than 5M® for an exterior 
perturber, 2Mq for an interior perturber in the 3:2 resonance, and < lM® for an interior perturber in the 2:1 resonance. 
These constraints are independent of the eccentricity of the perturber. For low initial eccentricity orbits (e ~ 0.01) the limits 
are as low as O.2M0, twice the mass of Mars. 

We marginalized over eccentricity to compute limits on the presence of coplanar planetary companions in a 2:1 and 1:2 
resonance with the transiting planet. Figure [7| shows the limits on planets within these resonances. Remarkably, the interior 
resonance can rule out at 3-<r the presence of earth-mass bodies of any eccentricity within 0.5% of exact resonance (we have 
not included bodies with eccentricities large enough that their Hill spheres overlap with that of HD 209458b). The exterior 
resonance is slightly less sensitive, but can still rule out ~ 2M® mass planets within 0.5% of resonance. We have checked 
for stability on timescales of several hundred orbits for resonant bodies, but have not carried out a longer timescale stability 
analysis. 
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4.1 Comparison with TrES-1 System 

The first analysis of this sort was conducted for the TrES-1 planetary system Jsteffen fc Ag O ]||2005h . That analysis used 11 
transits observed with ground-based telescopes and was also able to probe for sub-earth mass planets. Those observations had 
an average timing uncertainty near 110 seconds-the most precise being 26 seconds. For this analysis of HD 209458b, being 
observed by HST, the average timing uncertainties are significantly smaller, with an average of 25 seconds (17 seconds for the 
STIS data). 

The factor of four improvement in timing precision corresponds to a comparable improvement over the constraints on 
secondary planets in the HD209458 system. Figure |H| compares the limit on the mass of a companion planet for the TrES-1 
system with that of the HD209458 system for perturbers with an eccentricity of 0.02. This figure shows that the constraint 
on the HD209458 system is similar in shape but markedly lower. 

Though the space-based data can give tighter constraints for systems like TrES-1 and HD209458, figure|H|demonstrates the 
capacity that ground-based telescopes have to probe for secondary planets. The transit observations of the TrES-1 system were 
conducted on relatively modest telescopes (lm) and were not optimized for obtaining the time of transit. Even then, the best 
of the ground based observations (obtained using the 1.2m Fred Whipple telescope ICharbonneau et alj l2005) have a timing 
precision that is essentially equal to the average uncertainty of the HST data. Thus, additional ground-based observations, 
with precision comparable to the best of those obtained for the TrES-1 system could readily probe for companions with masses 
less than that of the Earth that are trapped in low-order, mean-motion resonances. 



4.2 Constraints Including Radial Velocity Measurements 

Theoretical studies l|Agol et ahfeoOsF) indicate that high-precision radial velocity measurements should provide more stringent 
constraints than TTV for planets that are not in resonance. We therefore conducted a combined analysis that utilize s both the 
transi t data and radial velocity data for this system. We used the radial velocity data that were published bv lLaughlin et all 
^0051) : discarding the data that they indicate were taken near transits. 

For this combined analysis we a pplied essentially the same algorithm that we used for the transit-only analysis and for 
the analysis of the TrES-1 system in lsteffen fc Agol (12005^ with the difference that we compare the radial velocity data to a 
model with two Keplerian orbits with the parameters of the two planets. The change in the orbital elements of the planets due 
to non-resonant planet-planet interactions are small enough to not affect the radial velocities, justifying the use of Kepler's 
equation. The results of our combined analysis ar e shown in Figure El We see from Fig ure El that our results are in good 
agreement with the theoretical predictions given in lAgol et alJ <l2005l) and lsteffen fc Agoll i2005l) . 



5 DISCUSSION 

This investigation started as part of a larger quest for detecting terrestrial-mass planets orbiting main-sequence stars. Transits 
times are particularly sensitive to earth-mass planets in mean-motion resonance with the transiting planet. For the case of the 
HD 209458 system, we can rule out a wide range of terrestrial-mass planets based on the HST data analyzed and presented 
here. The paradigm o f core accretion plus migration pre dicts that smaller planets ought to be trapped in resonance with larger, 
short-period planets Jzhou et aJbOQSt lThommesll2005l) . With the analysis of only two systems (HD 209458 and TrES-1), we 
cannot draw general conclusions, but in the next few years there will be likely be transit timing observations of all currently 
known transiting planets, so a significant constraint on the presence of resonant, earth-mass planets could be placed with this 
sample. In the process we have derived accurate parameters for the transiting planet which are broadly consistent with the 
analyses of other authors. 

Despite its lofty position above the Earth's atmosphere, it is apparent from these data that the Hubble Space Telescope 
is not ideally suited to timing transits. The main problem is its low-Earth orbit which causes occultation of the Sun leading 
to thermal variations within the orbit and causes occultation of the target star leading to gaps and re-pointing. The thermal 
variations cause sensitivity variations, while the gaps make a model for limb-darkening necessary to measure the transit times. 
If the data errors behaved according the Poisson statistics, the expected timing uncertainty for the STIS data would be of 
order 3-5 seconds, while the derived times are 4-5 times larger. The FGS data had a decreased sensitivity relative to the 
STIS data due to a lower count rate and the observation of each transit for only two or three orbits, an observing strategy 
which conserved Hubble time but was unsuited to the precise measurement of transit times. Future, precise measurements 
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of transit times are possibly better carried out with satellites in Earth-trailing orbit, such as the Spitzer Space Telescope. 
Spitzer has the extra advantage that the infrared is less subject to limb-darkening; however, it would collect fewer photons as 
it has a smaller aperture and is observing past the peak of the photon count spectrum of G-type stars. In addition, ground- 
based observations with accurate relative photometry and coverage from observatories separated in longitude should make 
very accurate transit timing feasible from the ground. Finally, the Kepler satellite will discover a large number of transiting 
jupiters - th e longer period planets w ill allow even more precise constraints on companions over a larger range of semi-major 
axis ratios jHolman fc Murray! Eooih . To contrain the presence of planets with mutual inclinations, more data points are 
required so that the model can accomodate more free parameters. 
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Figure 5. Constraints on the mass of an interior, secondary planet in the HD209458 system as a function of semi-major axis ratio and 
eccentricity of the secondary planet. We assume that the known planet has an initial eccentricity of zero. The contours correspond to 
100 (dotted), 10 (dashed), and 1 (solid) earth- masses. The dark region in the upper-left portion of the graph is where the orbits overlap. 
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Figure 6. Constraints on the mass of an exterior, secondary planet in the HD209458 system as a function of semi-major axis ratio and 
eccentricity of the secondary planet. We assume that the known planet has an initial eccentricity of zero. The contours correspond to 
100 (dotted), 10 (dashed), and 1 (solid) earth- masses. The dark region in the upper-left portion of the graph is where the orbits overlap. 
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Figure 7. Upper limit on the mass of any coplanar planetary companion present in the interior resonance (upper panel) or exterior 
(lower panel) 2:1 resonance with HD 209458b. The horizontal dashed line corresponds to an earth-mass plant. 
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Figure 8. A comparison of the limits on a secondary planets obtained by TTV analysis for the TrES-1 system with 11 transits (upper 
curve) and the HD209458 system with 13 transits (lower curve). The primary distinction between the transit data for each system is that 
the ground-based TrES-1 observations have a timing uncertainty that is a factor of four, on average, larger than those of the HD209458 
system, which were observed with HST. 
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Figure 9. Limit on the mass of a low-eccentricity secondary planet in the HD209458 system. The solid, black curve is obtained from an 
analysis of transit times and radial velocity measurements simultaneously assuming that the orbits are initially circular. The gray curve 
is the result of an analysis of the transit times only, with the perturbing planet initially having an initial eccentricity of 0.02. The thin 
dashed curve is from a perturbation theory calculation and are given by equations (A7) and (A8) in Agol et al. (2005). The large dots 
are from equation (33) in Agol et al. (2005), although without the factor of 4.5 which is the difference between the peak TTV amplitude 
and the TTV amplitude for intitially exact j:j + 1 resonance. The large, dashed line is from equation (2) in Steffen & Agol (2005) and 
is the expected limit from radial velocity measurements. The small, horizontal dashed line corresponds to the mass of the Earth. 
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